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Abstract 

A very intuitive description of nucleus-nucleus collision phenomena 
is provided by the relativistic fluid dynamics. We consider a l+f 
dimensional relativistic imperfect fluid flow to approximate the high 
energy heavy ion collision. The article investigates the application 
of the continuous symmetry group on the relativistic fluid energy- 
momentum tensor conservation equations in the ultrarelativistic limit 
7 — > oo. 

PACS: 25.75.-q, 47.75.+f 



1 Introduction 



There are three important theoretical approaches of the heavy ion collisions: 
fluid dynamics, Boltzman equation and statistical models. 

A simple but intuitive description of the general nuclear phenomena is 
given by fluid dynamics and it has been successfully applied to the ultrarel- 
ativistic central heavy ion collisions 

The Boltzman equation has already been applied to heavy ion collisions 
for relative low energies and there is also a relativistic Boltzmann equation 
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using Jiittner distribution y-'"**"'^^ = 1/ (27r/i)^ exp ((/i — p°'Mo)/T) where /i is 
the chemical potential, Ua and p° are the four- velocity and four-momentum 
of the particle, T is the temperature parameter 

In statistical models, one use the equilibrium thermodynamics to get 
the properties of the system. One of the main interest is the study of the 
phase transition from hadrons to quarks and gluons. At this point we can 
mention that the lattice quantum chromodynamics and subsequent use of 
Monte Carlo techniques enabled us to study physical observables over the 
entire temperature range from to oo |0] and the very existence of a phase 
transition. 

A consistent way of investigation of the high energy heavy ion collisions is 
the manifest covariant formulation of relativistic fluid dynamics. An excellent 
review about the status of fluid dynamics models for relativistic heavy ion 
reactions is the article written by D.Strottman One of the arguments to 
justify relativistic fluid dynamics is that the mean free path A is smaller than 
the size L of the nucleus. A simple estimate of the order of magnitude of A 
is about 1 fm, which is to be compared to the size of heavy nucleus 15 fm 
(for uranium). On the other hand, it was showed that a single particle 
quantum mechanics is formally analogous to the Euler hydrodynamics. 

Solving the hydrodynamics equations is hardly discouraged by the puzzle 
of choosing the rest frame. This situation leads us to acausal and instable 
solutions 0. Furthermore, the choose of the initial conditions is rather un- 
clear. In the same time the numerical approaches do not give us satisfactory 
quantitative results. 

Complementary information may be achieved by exploiting the Lie sym- 
metry group of the covariant relativistic hydrodynamics equations. The 
power of this technique consists in the possibility to explore the properties 
of physical systems, like the symmetry structure and the invariants, without 
solving the corresponding differential equations. The group structure and the 
invariants can help us to reduce the order of the equation and even to inte- 
grate them. Complex systems were successfully studied using this approach 
[0]. In we study a simple form of energy-momentum tensor conservation, 
where we already obtained symmetries and invariants. 
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2 Energy- momentum tensor conservation equa- 
tion 



Relativistic fluid dynamics is well describe by the number of particles N and 
energy-momentum tensor T^p conservation equations In the ideal case 
we have : 

Tap = PVaP + {P + P)UaUp 
Na = nUa 

where p is the pressure, p is the energy density, n is the number of particles 
density and Ua '■ (7/?, 7) is the 4-velocity field. 

There are two ways of choosing the rest frame : in the Landau way Ua 
is the energy transport velocity where Tjo = in the rest frame, while in the 
Eckart way Ua is the particle transport velocity where A^j = in the rest 
frame. The dissipation contribution is introduced by redefining the energy- 
momentum and number of particle tensor by adding correction terms : 

Tap = PVaP + {P + p)UaUp + AT^^/j, 

Na = nUa + AiV„ 

In the Eckart frame AA^^^ = 0, so the dissipation contribution is present 
only in the energy-momentum terms. In the following we choose the Eckart 
approach. The construction of the most general dissipation term /STap is 
coming up from the positivity of the entropy production: 

AT"^ = -'qH^'^H'^^W^s - x [h^'^U'^ - H^^U'') - (H'^^d^U"' 

where we have shear tensor: 

2 

Wap = dpUa + daUfi - -1]ap B^W , 

heat-fiow vector: 

Qa = daT + TU''dpUa, 

T is the temperature and projection tensor on the hyperplane normal to Ua 

HaP = VaP + UaUi3 . 
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We may identify x, C the coefficients of heat conduction, shear viscosity 
and bulk viscosity. 

Making some calculations we can write the energy-momentum tensor in 
the following form : 

Tap = Prial3+ {P + p)UaUi3 



daUp + dpUa - -iva/3 + Ua Up) d^W + Wd^{Ua Up) 



-x[UadpT+Upd^T + 2Ua UpWd^T + TWd^iUa Up)] 
-C {r]ap + UaUp)d^W 
With this form we rewrite the energy-momentum conservation : 



^df^iU^Up 



P + p+{-V-C]d-yU^-2xUW^T 



-a" [{ri + xT) Wd^{U^ Up) - ri {dJJp + dpU^) -x{UadpT+ Up T)] 

We can see this equation as a polinom in U, with the power 3,2,1,0 
. From this expression we select only the highest power of U and the free 
terms. Because U is proportional with 7 Lorentz U^ ^ U^ ^ U ^ 1 in 
the ultrarelativistic limit, or we can make a boost to a system were 7 is big 
enough to obtain to above inequality. This approximation is convenable for 
us because the form of the equation is much more simpler. In the same time 
we don't have any physical contradictions - as we will see later on Lorentz 
transformation will be one of the symmetry group transformation - and we 
don't loose physical information about the fluid. 

Doing this one can find the following form of the equation : 

d^T^p = dp\^p+ (^^Tj - C) d'^uj^ -2x-d^ {U^ UpUW^T) 

In the final form we have in the 1-1-1 dimensional approximation two 
equations: 

Px + {Iv - C) (uxx - v^t) - 2x {u^T^^ + uv'^Ttt - 2v?vT^t) = 
Pt + {In - C) {u.t - vtt) - 2x {u'vT,. + v^Tu - 2uv^T,t) = 



where Px,Pt, etc. means the partial derivative of p with respect to x,t, etc. 
The u and v are the spatial and temporal components of the velocity field 

U {'jv, 'j) = U {u, v). 



3 Symmetry group of differential equations 

The symmetry group of a system of differential equations is the largest local 
group of transformation acting on the independent and dependent variables 
of the system with the property that it transform solutions of the system to 
other solutions. 

We restrict our attention to local Lie group of symmetries, leaving aside 
problems involving discrete symmetries such as reflections. 

Let iS be a system of differential equations. A symmetry-group of the 
system 5 is a local group of transformations G acting on an open subset 
M of the space of independent and dependent variables for the system with 
the property that whenever u=f(x) is a solution of S, and whenever g ■ f is 
defined for g & G, then u = g ■ f{x) is also a solution of the system. 

Applying the standard procedure we solve the defining equations for 
the symmetry group of the given system of differential equations. We find 
a nine parameter symmetry group for the ultrarelativistic fluid dynamics 
equations. The basis of the corresponding solvable Lie algebra of this group 
is : 

Vi = dx (spatial translation) 

V2 = dt (temporal translation) 

V3 = Ot (temperature translation) 

V4 = dp (presure translation) 

V5 = xdr 
Ve = tdr 

Vj = udu + V dy — 2T Ot — pdp (dilatations) 
Vg = t dx + X dt — u dy — V du (hiperbolic rotations) 

= Ax dx + Atdt — udu — V dy + 2T St + ^p dp (dilatations) 
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The symmetry group infinitesimal generator is defined by : 
V = ^d, + Tdt + $9„ + ^d, + Tdr + ndp 

where 





= Ci + Cs ■ t + Acg ■ 


X 


r 


= C2 + Cs - X + Acg 


■ t 




= (C7 - Cg)^ - Cgf 






= (C7 - Cg)v - CsU 




r 


= C3 + C5X + C&t - 


2(C7 - Cg)T 




= C4 - (C7 - hCg)p 





and Oj are arbitrary group parameters. 

We calculate the second order prolongation of V : 

The coefficient functions of the prolongation of pr*^"^V are given by the fol- 
lowing formula : 

i=l 

where p is the number of the independent variables, ^* are the coefficients of 
the partial derivative of the independent variables (a;,t) - here ^* = (^,t) , q 
is the number of dependent variables u = {u,v,T,p) - in this case we have 
g = 4 so $Q = ($, r, Q) and Di is the total derivative 

df df 
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U 



J.i 



(9x* dx^dx^^...dx^k 
J = (ji, jjt), < #J < n , n is the highest order derivative appearing in 
/. For example <l>^ = D^i^ - ^u^ - rut) + iu^x + ru^t , ^"^"^ = Dl{m - ^v^ - 
TVt) + ^Vxxx + TVxxt and D^F = + r„M^. + T^v^ + T^T^ + Tpp^- 

An n-th order differential invariant of a group G is a smooth function 
depending on the independent and dependent variables and their derivatives, 
invariant on the action of the corresponding n-th prolongation of G []^ . 

We use the method of characteristics for computing the invariants for the 
prolongation of the every element of the Lie algebra ; then we re-express the 
next vector in terms of first vector's invariants; we repet this procedure until 
we obtaine all global invariants of the vectors. In the end we obtained seven 
global independent invariants of the symmetry group: 



h 
h 
h 



{u + v) {px + pt) exp 



2vx 



(m2 -t>2) (pl-pf) 



2\ V.t — Vx 



Ux - Vt 



PI - Pt 



{Px+Pt) {Ut - Vx) 



Pxx ~ Ptt 
2 I , \2 



{Pxx - Ptt) {U + V) 
{Px+Ptf {UxVt - UtVx) 

{UxVt - UtVx) {Px +Ptf 
{u + vf {Txx - Ttt) 



exp 



2v, 



Ux + Vt 



{UxVt - UtVxf {Txx - Ttt)' 



(2 1 2 

\Pxt ~ 2Pxx ~ PxxPtt 



[Pxt 2Pxx 



PxxPtt 



Uxt + UxtV X 



Vtxx {Utt + Vxx) 



( rp2 17^2 rp rp 

\-'-xt 2 XX -l-xx-l-tt 



Vlt + VxtUxx - \vl^ - Vxx {Uxt + Vtt) 



ult + UxtVxx 



1 2 

2^xx 



Vixx {Utt + Vxx J 
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4 Conclusions 



Relativistic imperfect fluid flow seems to be a very good approach of the 
ultrarclativistic heavy ion colhsion, because we can very weU suppose a zero 
mean free path and a instantaneous local equilibrium. 

In spite of the simplification of the equations, our results indicate a non- 
trivial structure of the symmetry group for the 1+1 dimensional relativistic 
fluid dynamics equations. 

As expected due to the covariant formulation of the theory, the symmetry 
group contain the Poincare generators. And of course, other system specific 
transformations mentioned above are present. 

Beside the physical meaning that can be associated to the invariants, they 
can be used for reducing the order of the original equations. Doing this one 
can hope to find simpler equations that can be integrated; unfortunately only 
for special type of groups one can find the general solutions of the equations 
by quadratures alone. 
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